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Abstract —In this paper, we consider decode-and-forward (DF) 
relay beamforming with imperfect channel state information 
(CSI), cooperative artificial noise (AN) injection, and finite- 
alphabet input in the presence of an user and J non-colluding 
eavesdroppers. The communication between the source and the 
user is aided by a multiple-input-multiple-output (MIMO) DF 
relay. We use the fact that a wiretap code consists of two parts: 
i) common message (non-secret), and ii) secret message. The 
source transmits two independent messages: i) common message 
(non-secret), and ii) secret message. The common message is 
transmitted at a fixed rate i?o, and it is intended for the nser. 
The secret message is also intended for the user but it should be 
kept secret from the J eavesdroppers. The source and the MIMO 
DF relay operate under individual power constraints. In order to 
improve the secrecy rate, the MIMO relay also injects artificial 
noise. The CSI on all the links are assumed to be imperfect 
and CSI errors are assumed to be norm bounded. In order to 
maximize the worst case secrecy rate, we maximize the worst 
case link information rate to the user subject to: i) the individual 
power constraints on the source and the MIMO relay, and ii) the 
best case link information rates to J eavesdroppers be less than 
or equal to i?o in order to support a fixed common message rate 
Ro. Numerical results showing the effect of perfect/imperfect CSI, 
presence/absence of AN with finite-alphabet input on the secrecy 
rate are presented. 

keywords: MIMO relay beamforming, physical layer security, multiple 
eavesdroppers, artificial noise, imperfect CSI, finite-alphabet input, semi- 
definite programming. 

I. Introduction 

The physically degraded discrete memoryless wiretap chan¬ 
nel model considered by Wyner in HI opened the path for 
reliable and secure information transmission using physical 
layer techniques. Subsequent extension to discrete memory¬ 
less broadcast channel and Gaussian channel was done in 
El and 0, respectively. A wireless network can be easily 
eavesdropped due to the broadcast nature of wireless transmis¬ 
sion. However, using physical layer techniques (e.g. wiretap 
codes, beamforming using multiple antennas, artificial noise 
injection etc.), a wireless network can be secured from getting 
eavesdropped. Achievable secrecy rate and capacity in single 
and multiple antenna wiretap channels have been reported by 
many authors, e.g., a-ii. 

A relay, operating in decode-and-forward (DE) or amplify- 
and-forward (AE) mode, can act as an intermediate node 
and help improving the secrecy rate HO). DE and AE relay 
beamforming techniques for secrecy under perfect/imperfect 
channel state information (CSI) have been well studied in the 
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literature, e.g. im-M- In these works, the transmit code¬ 
word symbols belong to an infinite constellation (Gaussian). 
However, in a practical communication system, the codeword 
symbols will belong to a finite alphabet set, e.g., M-ary 
alphabets. The effect of finite constellation on secrecy rate has 
been reported in lfT7l - E3l . In EH, DE relay beamforming for 
secrecy with finite alphabet has been considered. There it was 
shown that when the source power and relay beamforming 
vector obtained for Gaussian alphabet, when used with finite 
alphabet, could lead to zero secrecy rate. A power control 
algorithm was suggested to alleviate the loss in secrecy rate. 
Motivated by the above works, in this paper, we consider 
secrecy rate in DE relay beamforming with finite-alphabet 
input using a MIMO relay. The considered system consists of 
a source node, a destination node, and multiple non-colluding 
eavesdroppers. A DE MIMO relay aids the communication 
between the source and destination. It is known that secrecy 
rate can be improved through the use of artificial noise (AN) 
injection 0, ini, El, ED, El, El- In this work, we 
allow the MIMO relay to inject AN in addition to relaying 
the information symbol from the source. Consequently, we 
solve for both the optimum source power, signal beamforming 
weights as well as the AN covariance matrix at the MIMO 
relay. Since the CSI will not be perfect in practice, we consider 
a norm-bounded CSI error model and investigate the effect 
of imperfect CSI on the secrecy rate. We use the fact that 
a wiretap code consists of two parts: i) common message 
(non-secret), and ii) secret message. The source transmits two 
independent messages: i) common message (non-secret), and 
ii) secret message. The common message is transmitted at a 
fixed rate Rq, and its intended for the destination node. The 
secret message is also intended for the destination node but 
it should be kept secret from J eavesdroppers. The source 
and the MIMO DE relay operate under individual power 
constraints. In order to maximize the worst case secrecy rate, 
we maximize the worst case link information rate to the user 
subject to: i) the individual power constraints on the source 
and the MIMO DE relay, and ii) the best case link information 
rates to J eavesdroppers be less than or equal to Rq in order to 
support a fixed common message rate Rq. Numerical results 
showing the effect of perfect/imperfect CSI, presence/absence 
of AN with finite-alphabet input on the secrecy rate are 
presented. 

Notations : A £ cNixATj jj^pHes that A is a complex 
matrix of dimension Ny x N2. A ^ 0 and A 0 imply 
that A is a positive semidefinite matrix and positive definite 
matrix, respectively. Identity matrix is denoted by I. Transpose 
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Fig. 1. System model for MIMO DF relaying. 

and complex conjugate transpose operations are denoted by 
and [.]*, respectively. E[.] denotes expectation operator. 
||.|| denotes 2-norm operator. Trace of matrix A G is 

denoted by Tr(A). G ^ CM(0, ’S') implies that is 

a circularly symmetric complex Gaussian random vector with 
mean vector 0 and covariance matrix ’S'. 

II. System model 

Consider a DF cooperative relaying scheme which consists 
of a source node S having single transmit antenna, a MIMO 
DF relay node i? having N receive/transmit antennas, a 
destination node D having single receive antenna, and J non¬ 
colluding eavesdropper nodes Ei,E 2 ^--- ,Ej having single 
receive antenna each. The system model is shown in Fig. [T] 
In addition to the links from relay to destination node and 
relay to eavesdropper nodes, we assume direct links from 
source to destination node and source to eavesdropper nodes. 
The complex channel gain vector between the source and the 
relay is denoted by g = [ 51 , 32 , ■ ■ ■ ,5Af]^ S Likewise, 

the channel gain vector between the relay and the destination 
node D is denoted by h. = [hi, / 12 , • • ■ , G and the 

channel gain vector between the relay and the 5 th eavesdropper 
node Ej, 1<J< J, is denoted by Zj = [zij, Z 2 j, - ■ ■ , S 
(j^ixiv rpjjg channel gains on the direct links from the source to 
D and the source to Ej are denoted by hg and zoj, respectively. 

The MIMO relay operates in half duplex mode, and 
the communication happens in two hops. Each hop is di¬ 
vided into n channel uses. We use the fact that a wire¬ 
tap code consists of two parts: i) common message (non¬ 
secret), and a) secret message. In the first hop of trans¬ 
mission, the source S transmits two independent messages 
Wo and Wi which are equiprobable over {1,2, ••• 
and {1,2,-•• ^ respectively. Wq is the common 

message which is transmitted at a fixed rate Rq and its 
intended for the destination D. Wi is a secret message which 
is transmitted at some rate Rs{Ro) and its also intended 
only for D and it should be kept secret from all EjS. For 
each Wq and Wi drawn independently and equiprobably 
from the sets {1,2,--- and {1,2,-•• 


respectively, the source S maps Wq and Wi to a codeword 
{xm}m=i of length n. Each symbol, Xm, in the codeword is 
independent and equiprobable over a complex finite-alphabet 
set A = { 01 , 02 ,-•• jQm} of size M with E[a:m] = 0, and 
E[|a:m|^] = 1. The source is constrained by the available power 
Ps and it transmits the weighted symbol which is y/P^Xm in 
the mth channel use, where 1 < m < n, and 0 < Ps < Ps- 
Hereafter, we will denote the symbol Xm of the codeword 
{xTn}m=i by 2 :, and we will consider only one channel use. 

Let y^, yjj, and yE^j denote the received signals at the 
MIMO relay R, destination D, and jth eavesdropper Ej, 
respectively, in the first hop. We have 


Vr 

= y/Kgx + J7r, 

(1) 

yoi 

= y/p'shox + rjDi, 

(2) 

VEii 

= \/~PsZ0jX -b PBi,, 

(3) 


where CM{^,NoI)), CA/'(0, Wq)), and 

CAf{0, Nq)) are receiver noise components and are assumed 
to be independent. 

In the second hop of transmission, MIMO relay applies 
the complex weight cj) = [(j>i, 4 > 2 , ■'' ^ on the 

successfully decoded symbol x and retransmits it. In order to 
improve the secrecy rate, MIMO relay also injects the artificial 
noise xp G CA/’(0,^)). The symbol transmitted by 

the MIMO relay on the ith, 1 < z < antenna is piX + ipt. 
Let 5 D 2 > yE 2 j denote the received signals at the destination 

D, and jth eavesdropper Ej, respectively, in the second hop. 
We have 

yD2 = hcpx + h-ip + riD2, ( 4 ) 

yE2j = ZjCpx -b Zjip -b r]E2E (5) 

where CAf{0,No)), and riE 2 j{^ CM{0,No)) are re¬ 

ceiver noise components and are assumed to be independent. 
Using (l2]i, (|4]i, and Q, Q, we rewrite the received signals at 
D and Ej in the following vector forms, respectively: 

Vd = [j/Di, yD 2 f 

'J-' 

= iV^ho, hep] X + [r]Di, hi> + r]D2f, (6) 
Vei = {yE,E yE2jf 

r— T 

= [y/WsZoj, Zj(p] X -b [yEii, Zjip + r]E2jfi^) 

We assume that the MIMO relay’s transmit power, denoted 
by Pr, is constrained by the available power Pr. This implies 
that 

Pr = E{\\{ePx + 

= \\cpf + Tr(’®') < Pr. (8) 

We also assume that the channel remains static over the 
entire codeword transmit duration. Eurther, denoting the secret 
message decoded at the MIMO relay R and destination D by 
WA and WA, respectively, the reliability constraints at R and 
D and the perfect secrecy constraints at P^s are as follows: 

¥AWAAWi) < e„, ¥AWA AWi) < e„, 

^I{Wi-y%) < e„, Vj = 1,2,--- ,J, 
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where is the received signal vector at Ej in 2n channel 
uses, and —>■ 0 as n —>■ oo. We also note that the reliability 
constraints at the MIMO relay R and destination D for the 
secret message also ensure the reliability of the common 
message. 


III. DF RELAY BEAMFORMING - PERFECT CSI 


In this section, we assume that the CSI on all the links are 
known perfectly. Using ([TJ, (|6]l, and (|7]i, we get the S-R, S-D, 
and S-Ej link information rates, respectively, as follows: 


where 


2 V ^0 y’ 

1 (Ps\K\^ hct)4>*h* \ 

2 \ No ^No + h^h*)' 

^j( Ps\zojf \ 

2 'v No No + z,'Sfz*J 


M 


= TtX! / Pn{y - '/po-l) 


log 


M 

1^1 

Pnijj - y/pai) 

2 M 

JJ S Pniy ~ y/pO-m) 

m—1 


-dy, 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


and PniP) = ■ The factor 1/2 in (|9l), ( fTOl i, and (fTTl i 

is due to two hops. Further, the MIMO relay R will be able 
to decode the symbol x if the following condition holds true 
0111151161221 : 

> ^I{x;yjo). (13) 

In order to find the maximum achievable secrecy rate Rs{Ro) 
which also supports the fixed common message rate Ro, we 
maximize the S — D link information rate subject to i) S — Ej, 
1 < j < J, link information rates be less than or equal to Ro, 
ii) the information rate constraint in (fT^ . and in) the power 
constraints. The optimization problem is as follows: 

Rd{Ro) = max l-I{x;yjj) (14) 

Ps, 2 

s.t. ^I{x;yE.) < Ro, Vj = 1,2, • • • , J, (15) 

^I{x-,yR) > ^I{x;yjo), (16) 

0<Ps<Ps, 0, + Tr(^) < Pr. (17) 


Having obtained Rr,{Ro) from (fT4l i. the maximum achievable 
secrecy rate Rs{Ro) for a given common message rate Rq is 

0 

RsiRo) = {RDiRo) - Ro}^, (18) 

where {a}^ = max(0,a). From the constraint in (fThl l. its 
obvious that the upper bound for S — D link information rate, 
denoted hy Rd, can be obtained by evaluating 0 at Pg = Ps- 


For the values of Ro over the interval [0, Rd], the maximum 
achievable secrecy rate, denoted by Rg, is obtained as follows: 

Rs = max {i?_D(i?o) - i?o}^ (19) 

0 < Ro < Rd 

= max {RrUAi) — (20) 

0 < i < L 


where L is a large positive integer, Ai = Rr,/L, I is an 
integer, and Ro = /Ai. 

We solve the optimization problem (fT4li for a fixed Pg = 
kA 2 , where A 2 = Pr/K, K is a large positive integer, 
and 1 < k < K. Hereafter, we will assume that Pg is 
known. Further, it is shown in 1271281 that for various M-ary 
alphabets, mutual information expression in (fTSli is a strictly 
increasing concave function in SNR. With this fact, we rewrite 
the optimization problem (fl4li into the following equivalent 
form: 


max 


/ h^h* \ 

No + h'^h* ) 


( 21 ) 



s.t. Vj = 1,2,--- , J, 


No + Zj^z*J 


< I-H 2 R 0 ), 


c > 


hi^h* 


No + h^h 

A 0, rank{^) = 1, (F A 0, Tr($ + ’!')< Pr, 


( 22 ) 

(23) 

(24) 


where = cjicj)*, a = (-^4^), bj = ), and 

c = Further, relaxing the rank{^) = 1 constraint, 

we rewrite the above optimization problem into the following 
form: 



max 

t 

(25) 


t, 

4” 


s.t. 

(t - a)(No + h'^h*) - (h^h*) 

< 

0, 

(26) 


Vj = 1,2,... 




(z,^z*) 

-(l-^(2Ro)-b,)(No + z,^z*) 

< 

0, 

(27) 


(h^h*) - (c - a)(No + h'^h*) 

< 

0, 

(28) 


$ A 0 , if A 0 , Tr(# + ’S') < Pr. 

(29) 


The above problem can be easily solved using bisection 
method 1291 . The initial search interval in the bisection method 
can be taken as [0, c]. In the appendix, we show that the 
solution $ of the above problem has rank 1. Further, denoting 
the maximum value of t by tmax, the secrecy rate is obtained 
as follows: 

RsiRo) = "I — .Roj’ • (30) 

IV. DF RELAY BEAMFORMING - IMPERFECT CSI 

In this section, we assume that each receiver has perfect 
knowledge of its CSI. We also assume that the control unit 
which computes the source power, signal beamforming vector 
and AN covariance matrix has imperfect CSI on all links. The 
imperfection in CSI is modeled as follows 1151251261 : 

g = g + Eg, ho = ho + Cho, h ~ h + er, 

Vj = 1,2, • • • , J, zoj = Zoj ^ZQj, A , (31) 
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where g, Hq, h, 'z^j, 'Zj are the available CSI estimates, and 
Bg, Cho, Gh, Bzgj, Bzj are the corresponding CSI errors. We 
assume that the CSI errors are bounded, i.e., 

||6g||< Cg, |e?io| < e/toi ||ch.||^ etii 

Vj = , J, \ezo^\<€zg^, \\ezj<ez^. (32) 

With the above CSI error model, we write the rank relaxed 
optimization problem (l2Tl i as follows: 

/ uNq + (/r + eh,)(a^ + €•) {h + eh)* 

max min -^- 

\ No + ih + eh)'^{h + eh)* 



s.t. ^ 0, ’i' h 0, Tr($ + ’4') < Pr, (43) 

Veh s.t. Ile^ll^ <4i^ 

0 < ri < aWo + (^ + e^)(a'i'+ ^)(^ + e^)*, (44) 

Ve^ s.t. Ile^ll^ < 4 ^ 

Nq +(Ji +eh)^{h +Bh)* < r2, (45) 

^ < r\2Ro), Vj = l,2,...,J (46) 

S2j 

s.t. \\ez^t <el. 

bjNo + {zj + ezj)(bj'i' + ^){zj + e^^)* < sij, (47) 
s.t. \\ezj^<el. 

0 < S 2 j < No + (Zj + Bzj)'^ {Zj + Bzj)*, (48) 


s.t. 

Vj = l,2,--. ,J. 

63 iVo + (2j+e^^.)(63 «'+■*•) 

< I-\2Ro), 


II2 ^ 2 

^h\\ < 


S.t. 


c > max ( am 

V 

s.t. 


— ’ 

A ^0 + (^+eh.)'4^(^+eh.)* 

Ilehll" < 4. 


$ ^ 0, ^^0, Tr($ + ^') < Pr, 


(34) 

(35) 

(36) 

(37) 


where 


a = 


' Ps\\ho\ - Choi 
No 


' Ps\ ||g|| -eg 
^ No 


if {\ho\>eho), 0 else, (38) 


b,= 


^sll^ojl + e. 


Zoj I 


i 

-), (39) 


-) if (Il 9 ll> £ 9)5 0 else, (40) 


(^max — 


-Ps||(to| + efip 

iVn 


-). (41) 


The objective function in (l3^ corresponds to the worst case 
S — D link information rate over the region of CSI error 
uncertainty. The constraint in dTSl l corresponds to the best 
case S — Ej link information rate over the region of CSI 
error uncertainty. The constraint in (l3^ is associated with 
the information rate constraint in (|2^ . i.e., the worst case 
information rate to the MIMO relay R over the region of CSI 
error uncertainty should be greater than or equal to the best 
case information rate to destination D. 

Solving the optimization problem (l3^ is hard due to the 
presence of br in both the numerator and denominator of 
the objective function in (l33l l and the constraint in (l36l l. 
Similarly, appears in both the numerator and denominator 
of the constraint in (l35l l. So, by independently constraining the 
various quadratic terms appearing in the objective function in 
(l3^ and the constraints in (l35T l. (l36] l. we get the following 
lower bound for the above optimization problem: 


f’l 


T"! , ^2 , r3 , T’ 4 , 

^Ij. '* 2 j> 4 = 1 , 2 ,--- 


r2 


( 42 ) 


0 ^ {^vaax ^ ^ ’ 

(49) 

Ve^ s.t. ||eh4 <(1.=^ 


(h + eh)^(h + eh)* < ro, 

(50) 

Ve^ s.t. ||eh4 <(1.=^ 


< No + (h + eh)^(h + eh)*. 

(51) 


The quadratic inequality constraints in (l44l i and (l45l l are 
associated with the objective function in (l33]) . The constraint 
in (l46l l, and the quadratic inequality constraints in ( l47b and 
(|48] | are associated with the constraint in dTSl ). Similarly, the 
constraint in ( |49] |. and the quadratic inequality constraints in 
(I 50 I 1 and (ISTT i are associated with the constraint in (l36l l. 

Further, using ^-procedure ll^ . we transform the quadratic 
inequality constraints in (l44l i. (l45l) . (l47l) . (l48l l. djol l. and 
dSD, into the following linear matrix inequality (LMI) forms, 
respectively: 


ri > 0 , Ai > 0 , Ai = 

(a^ + ^) + Ai J (a^ + 

h(^a'^ + aNo + h.(a^ + — ri — Aie^ 

A2 > 0, A2 — 

[ —h^* —Nq — h^Sfh* + r2 — A2e^ J 

+ ^)'' -bjNo+‘^)z* + sij - 

^2j ^ 0 , l-L2j ^ 0 , ^2j — 

^ + tl2jl '^Z* 

Zj^* Nq — S 2 j — 

A3 > 0, A3 — 

r -^ 4 -A 37 -^h* 1 

L -h^* -h^h* J’ 

r4 >0, A4 > 0, A4 — 

r ^ + a4/ ^h* 1 

!_ Nq- r 4 - X 4 el^ J’ 

where Ai >: 0 , A 2 >: 0 , A 3 ^ 0 , A 4 ^ 0 , Bij ^ 0 , 

B2j y 0. We substitute the above LMI constraints in the 
optimization problem (|42]) . We get the following equivalent 
form for the optimization problem (|42]) : 
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max r (52) 

«i», 'I' I 

r’i, - r 4 , Xi,--- ,X 4 , 
sij, S2j, Mlj, M2j. j = l,2,- -,J, 

s.t. ^^0, ^^0, Tr($ + ^') < Pr, 
rr2 - ri < 0, Vj = 1, 2, • • • , J, sij - S2jl~^{2Ro) < 0, 
n >0, r4> 0, Ai >0, A 2 > 0, As >0, A 4 > 0, 

Ai to, A2t 0, As ^0, 244 ^ 0, 

S2j > 0, fiij > 0, fi2j > 0, Bij t 0, B2j t 0, 

rs- {c- amax)r4 < 0. (53) 

The above problem can be solved using the bisection method 
as discussed in Section mil The initial search interval in the 
bisection method can be taken as [0, c], where c is as defined 
in (Eil. Further, denoting the maximum value of r by Vmax, 
the lower bound on the secrecy rate is obtained as follows: 

Rs{Ro) > |i/(w)-i?o| • (54) 

V. Results and discussions 

In this section, we present numerical results on the secrecy 
rate for BPSK alphabet (i.e., M = 2), with/without AN, 
perfect/imperfect CSI conditions. We assume that N = 2, 
J = 1,2, 3, No = l, Ps = 0 dB, and Pr = 9 dB. 

Perfect CSI case of Section\nI\: We have used the following 
channel gains in the simulations: 


g = 

[-0.5839 -f 2.2907i, -0.7158 -f 0.1144i]^ 

,(55) 

ho = 

-0.3822- 0.3976i, 

(56) 

Zqi = 

0.0123-f0.0137z, 

(57) 

Z02 = 

0.0231- 0.0178Z, 

(58) 

^03 = 

-0.0045 - 0.0042i, 

(59) 

h = 

[0.2174 - 0.6913i, -0.4047 - 0.3159i[, 

(60) 

Zi = 

[0.3826 -f 0.0811i, 0.8389 - 0.0943i], 

(61) 

Z2 = 

[0.2977 -f 0.7902i, -0.2069 -f 0.4696i[, 

(62) 

Z3 = 

[-0.6076 -f 0.6637Z, -0.3316 -f 0.192H]. 

(63) 


In Fig.|2] we plot the secrecy rate versus Rq for BPSK alphabet 
(i.e., M = 2), with/without AN, J = 1,2,3 eavesdroppers, 
Ps =9 dB, and Pr = 9 dB. We observe that the secrecy rate 
initially increases with increase in Rq and then drops to zero 
for large values of Rq. We also observe that the injection of 
AN improves the secrecy rate when J = 2, 3 eavesdroppers 
are present. However, when only one eavesdropper is present, 
the secrecy rate plots with/without AN overlap. This is due 
to the null signal beamforming by the MIMO relay at the 
eavesdropper. This is possible only when the number of 
eavesdroppers is strictly less than the number of antennas 
in the MIMO relay which happens to be true for this case 
with N = 2 and J = 1. For the case when J = 1, the 
secrecy rate maximum happens at Rq = 0.001445. Further, 
for the case when J = 2,3 and without AN, the secrecy 
rate maximum happens at Rq = 0.145797, and with AN it 
happens at Rq = 0.080959 and Rq = 0.099059, respectively. 
It is seen that the secrecy rate falls approximately linearly for 



Fig. 2. Secrecy rate vs Rq in MIMO DF relay beamforming for BPSK 
alphabet, with/without AN signal. N = 2, Nq = 1, J = 1,2,3, M = 2, 
fixed Ps = 0 dB, and Pji = 9 dB. 

large Rq. The near-linear fall in secrecy rate for large values 
of Rq is due to the saturation of S — D link information rate to 
^ log 2 2 = 0.5 for M = 2. We have also numerically observed 
that the rank of $ is 1. 

Imperfect CSI case of Section |73 .• Here, we assume that 
the channel gains in (l55]) - (l63]) are the available CSI estimates. 
We also assume that the magnitudes of the CSI errors in all 
the links are equal, i.e., Cg = Chg = Czoj = = e. 

We solve the optimization problem (l52l i for BPSK alphabet 
(i.e., M = 2), with AN, fixed Rq = 0 . 0810 , Ps = 9 dB, and 
Pr = 9 dB. In Fig. [3 we plot Rs vs e with AN for J = 1, 2,3. 
We observe that the secrecy rate decreases with increase in CSI 
error and with increase in number of eavesdroppers. We have 
also numerically observed that the rank of $ is 1. 

VI. Conclusions 

We considered MIMO DF relay beamforming with im¬ 
perfect CSI, cooperative artificial noise injection, and finite- 
alphabet input in the presence of an user and multiple non¬ 
colluding eavesdroppers. The source transmits common and 
secret messages which are intended for the user. The common 
message is transmitted at a fixed rate Rq. In order to maximize 
the worst case secrecy rate, we maximized the worst case link 
information rate to the user subject to: i) the individual power 
constraints on the source and the MIMO DF relay, and ii) 
the best case link information rates to J eavesdroppers be 
less than or equal to Rq in order to support a fixed common 
message rate Rq. Numerical results showing the effect of 
perfect/imperfect CSI, presence/absence of AN with finite- 
alphabet input on the secrecy rate were presented. We would 
like to remark that the work presented in this paper can be 
extended to amplify-and-forward relay channel. 
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Fig. 3. Rs vs e in MIMO DF relay beamforming for BPSK alphabet and 
with AN signal. TV = 2, TVq = 1, = 1, 2, 3, M = 2, fixed Rq = 0.0810, 

Ps =0 dB, and = 9 dB. 

Appendix 

In this appendix, we analyze the rank of the optimal solution 
$ obtained from the optimization problem ( |25] ). We take the 
Lagrangian of the objective function —t with constraints in 
(I26ll-(l29ll as follows EH: 

4>, 'S', A, Ai, A 2 , A*, 0 = -f-Tr(Ai$) 

- Tr(A 2 'i') + a( Tr($) + Tr('5') - Pn) 

- a)(Ao + h^h*) - (h^h*)'^ 

J 

+ - (/-'(2i?o) - b,) {No + 

i=i 

+^(^{h^h*) - {c- a){No + (64) 

where A > 0, Ai ^ 0, A 2 ^ 0, p. > 0, Vj > 0, C ^ 0 
Lagrangian multipliers. The KKT conditions for ( l64l l are as 
follows; 

(al) all constraints in (l26])-(l29b, 

(a2) Tr(Ai$) = 0. Since Ai ^ 0 and ^ 0 Ai4> = 

0 , 

(a3) Tr(A 2 'i') = 0. Since A 2 ^ 0 and 'i' ^ 0 A 2 ’®' = 

0 , 

(a4) A(Tr($) + Tr(’i')-Pfl) = 0, 

(a5) fi(^{t-a){No + h^h*)-{h^h*)'^ = 0, 

(a6) Vj = l,2,... ,J, z/,((^,$^*)-(/-i(2i?o)-6,)(iVo+ 

z,^z*)) = 0 , 


(a7) S,(^{h^h*)-{c-a){No + h^h*)'j = 0, 

(a8) = 0 => p(Ao + = 1. This further implies 

that p > 0, 

(a9) §i = 0 Al = XI - fj.{h*h) + J2j=i’^jiZjZj) + 

ah*h), 

(alO) ^ = 0 A 2 = AT + p(f - a)ih*h) - 

- bj^ {z*Zj) - ^(c - a)(h*h). 

The KKT conditions (a8) and (a5) imply that the constraint 
( l26b will be satisfied with equality. Assuming $ 0, this 

further implies that t > a. The constraints dZTl i and ( l28l l imply 
that I~^{2Ro) > bj and c > a. The KKT conditions (a9), 
(alO), (a2), (a3), (a4), (a5), (a6), and (a7) imply that 

j 

XPr - ti{t - a)No + ^ (/"^(2i7o) - b.'^No 

+C(c - a) Ao = 0. (65) 

Let Pr be small enough such that the constraint in (|28]) is 
satisfied with strict inequality. This implies that the KKT 
condition (a7) will be satisfied only when ^ = 0. With •J = 0. 
we consider the scenario when the expression ( |65] | is satisfied 
for A > 0. With A > 0, the KKT condition (a4) implies that 
Tr($) + Tr(’®') = Pr, i.e., the entire relay power, Pr, will be 
used for transmission. Further, we rewrite the KKT condition 
(a9) as follows; 

J 

Al + fj,{h*h) = XI + ^ Vj{z*Zj) >- 0. (66) 

i=i 

The above expression implies that rank{Ai) > N — 
rank{^{h*h)) = N — 1. The KKT condition (a2) further 
implies that rank{Ai) = A — 1 and rank{^) = 1. 

We now show that the solution of the optimization 
problem ( l25l l has rank-1 even for large values of Pr. Let 
$ ^ 0 0) and ^ ^ 0 0) be the optimal solutions 

of (I 25 I 1 with 

Tr($) + Tr(’S') = P < Pr, 

and the objective function value t > 0. Dehne 

$ $ 

° “ Tr($) + Tr(’®') “ P’ 

tF ’3> 

° ^ Tr($) + Tr(’®') ^ P’ 

It is obvious that the objective function value, t, in the opti¬ 

mization problem (|25]) is a non-decreasing function in Pr. As 
discussed previously for small values of Pr, the optimization 
problem (|25] | attains it’s maximum value when entire power is 
used, i.e., ($, '5') == (P^o, -f^o) = {Pr^o, Pr^o)- This 
implies that the objective function value, t, in (|25]) is a strictly 
increasing function in Pr for small values of Pr. We now 
fix the directional matrices ($oi ’®'o) which are obtained for 
small values of Pr such that the constraint in ( |28] ) is satisfied 
with strict inequality. We rewrite the constraints in dZTl ) and 
(|28]) in the following forms, respectively; 
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Vj = ,J, r\ 2 Ro) > 


Zj^Z* 


No+zj'^z* 


c > { a + 


h^h* 


No + h^h* 


(67) 

( 68 ) 


In the above inequalities, the derivatives of the functions 
( jVo+y^.» ) w.r.t. P when evaluated at 

(P$o, P^o) are > 0 and > 0, respectively. This implies 
that the right hand sides of the inequalities in dhTl) and (l68l l 
are non-decreasing and strictly increasing functions in P, 
respectively, at (P$o, Pd^o)- This further implies that if 
the above inequalities are satisfied at {Ph^q, Pji^o), the 
optimization problem dZST l will attain its maximum value at 
(Pr#0: Pr'^o), i-e., when the entire available relay power, 
Pr, is used. When Pr is large such that the above inequalities 
fail to satisfy at {Ph^q, Pji^o), the optimization problem 
(IZSl l will attain its maximum value at (P^O; P^o)^ where 
P (< Pr) is the maximum power at which the above 
inequalities are satisfied at (P$0: The excess power 

(Pr — P) will remain unused. This implies that the ranks of 
and remain constant for large values of Pr. 
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